Lepton flavor violating Higgs boson decays 
in the MSSM-seesaw 
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Lepton flavor violating Higgs boson decays (LFVHD) are studied in the context of the Minimal 
Supersymmetric Standard Model (MSSM) enlarged with three right handed neutrinos and 
their supersymmetric partners, and with the neutrino masses being generated by the seesaw 
mechanism. We compute the partial widths for these decays to one- loop order and analyze 
numerically the corresponding branching ratios in terms of the MSSM and seesaw parameters. 
We analyze in parallel the lepton flavor changing lj — > U<y decays and explore the maximum 
predicted rates for LFVHD, mainly for H°,A a — > tJx decays, by requiring compatibility with 
neutrino and BR(lj — » U'y) data. We find LFVHD ratios of up to 10~ 5 in some regions of the 
MSSM-seesaw parameter space. 



1 Introduction 

The observed neutrino masses do require a theoretical framework beyond the Standard Model 
of Particle Physics with just three massless left-handed neutrinos. Within the MSSM-seesaw 
context, which will be adopted here, the MSSM particle content is enlarged by three right 
handed neutrinos plus their corresponding supersymmetric (SUSY) partners, and the neutrino 
masses are generated by the seesaw mechanism. Three of the six resulting Majorana neutrinos 
have light masses, m Ui ,i = 1, 2, 3, and the other three have heavy masses, m^, i = 1, 2, 3. These 
physical masses are related to the Dirac mass matrix m £> , the right-handed neutrino mass matrix 
iriM, and the unitary matrix Umns by diag{m Vll m U2 ,m U;i ) ~ Uj {NS {—mDm^jrn^ ) )UMNS arm 
diag{mN 1 ,mN 2 -,mN 2 ) ~ tum, respectively. Here we have chosen an electroweak eigenstate basis 
where mu and the charged lepton mass matrix are flavor diagonal, and we have assumed that 
all elements in mo = Y v < H% >, where Y v is the neutrino Yukawa coupling matrix and 
< i?2 >= wsin/? (y = 17A GeV), are much smaller than those of tum- The two previous 
relations can be rewritten together in a more convenient form for the work presented here as, 
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UT IN o, where R is a general complex and orthogonal 3x3 matrix, 



which will be parameterized by three complex angles 0j, i = 1,2,3. 

One of the most interesting features of the MSSM-seesaw model is the associated rich phe- 
nomenology due to the occurrence of lepton flavor violating (LFV) processes. Whereas in the 
standard (non-SUSY) seesaw models the ratios of LFV processes are small due to the small- 
ness of the light neutrino masses, in the SUSY-seesaw models these can be large due to an 
important additional source of lepton flavor mixing in the soft-SUSY-breaking terms. Even 
in the scenarios with universal soft-SUSY-breaking parameters at the large energy scale asso- 
ciated to the SUSY breaking Mx, the running from this scale down to tum induces, via the 



neutrino Yukawa couplings, large lepton flavor mixing in the slepton soft masses, and provides 
the so-called slepton- lepton misalignment, which in turn generates non-diagonal lepton flavor 
interactions. These interactions can induce sizable ratios in several LFV processes with SM 
charged leptons in the external legs, which are actually being tested experimentally with high 
precision and therefore provide a very interesting window to look for indirect SUSY signals. 
They can also induce important contributions to other LFV processes that could be meassured 
in the next generation colliders, as it is the case of the MSSM Higgs boson decays into r/2, re 
and /ue which are the subject of our interest. 

Here we compute the partial widths for these lepton flavor violating Higgs boson decays 
(LFVHD) to one-loop order and analyze numerically the corresponding branching ratios in 
terms of the MSSM and seesaw parameters, namely, Mq, M1/2, tan/3, m^y. and R. We analyze 
in parallel the lepton flavor changing L — > 1^ (i ^ j) decays and explore the maximum predicted 
rates for LFVHD, mainly for H , A — > r/2 decays, by requiring compatibility with BR(lj — > l^) 
data. For these we use the present experimental upper bounds given by ^ \BR(fi — » cy) | < 
1.2 x HT 11 , \BR{t -> pi)\ < 3.1 x 1(T 7 and \BR(r -> &y)\ < 2.7 x 1(T 6 . 

For the numerical analysis we choose, Mx = 2 x 10 16 GeV and Aq = 0. The Umns 
matrix elements and the m Vi are fixed to the most favored values by neutrino data^ with 

^Am^ = 0.008 eV, ^Am 2 atm = 0.05 eV, 9 12 = sol = 30°, 9 23 = Q atm = 45°, 13 = 0° and 
= a = (3 = 0. We consider two plaussible scenarios, one with quasi-degenerate light and 

A Am.^ 

degenerate heavy neutrinos and with m Ul =0.2 eV , m V2 = m Ul + 2m soi , m U3 = m Ul H — 2m atm 
and mjVi = mj\r 2 = ra^ 3 = rriN] and the other one with hierarchical light and hierarchical heavy 
neutrinos, and with m vi ~ OeV ,m U2 = J Am^ , m V3 = \J Am^ m and < niN 2 < rriN 3 - 

This is a reduced version of our more complete work'"' to which we address the reader for 
more details. 



2 Numerical results and conclusions 



We show in figs. Q through Q) the numerical results for the branching ratios of the LFVHD 
together with the branching ratios for the relevant lj — > lij decays. The results of BR{Hq — > r/2) 
as a function of tun, for degenerate heavy neutrinos and real R, are illustrated in fig. (^1), for 
several tan/3 values, tan /3 = 3, 10, 30, 50. Notice that in this case, the rates do not depend on R. 
The explored range in is from 10 8 GeV up to 10 14 GeV which is favorable for baryogenesis. 
We also show in this figure, the corresponding predicted rates for the most relevant lepton decay, 
which in this case is fj, — > ey, and include its upper experimental bound. We have checked that 
the other lepton decay channels are well within their experimental allowed range. The ratios for 
Aq decays, not shown here for brevity, are very similar to those for Hq decays in all the studied 
scenarios in this work. We have also found that the ratios for the light Higgs boson, ho, behave 
very similarly with ttin and tan /3 but are smaller than the heavy Higgs ones in about two orders 
of magnitude. ;From our results we learn about the high sensitivity to tan/3 of the LFVHD 
rates for all Higgs bosons which, at large tan /?, scale roughly as (tan /?) 4 , in comparison with the 
lepton decay rates which scale as (tan /3) 2 . The dependence of both rates on is that expected 
from the mass insertion approximation, where BR(H X — ► Ijli), BR(lj — > ^7) oc \mj^ log(mAr)| 2 . 
We find that the largest ratios, which are for Hq and Aq, are in any case very small, at most 
10~ 10 in the region of high tan/3 and high mjv- Besides, the rates for \i — > ey decays are below 
the upper experimental bound for all explored tan/3 and values. The branching ratios for 
the Higgs boson decays into re and /xe are much smaller than the r/2 ones, as expected, and 
we do not show plots for them. For instance, for = 10 14 GeV, and tan /3 = 50 we find 
BR(HW -> rp)/BR(Hto -> re) = 4 x 10 3 and BR(H^ -> rp,)/BR(H^ -> /ie) = 1.2 x 10 6 
for the three Higgs bosons. 



BR(H -> r/2) BR(n 67) 




Figure 1: Dependence of LFV ratios with ttiat (GeV) for degenerate heavy neutrinos and real R and for several 
values of tan/3, (la) BR(Ho — * Tj£). (lb) BR(n — > e7). The horizontal line is the experimental upper bound. 
In both plots, the solid, dashed, dashed-dotted and dotted lines are the preditions for tan/3 = 3, 10,30 and 50, 

respectively, and Mo = 400 GeV, M 1/2 = 300 GeV. 




Figure 2: LFV ratios for hierarchical heavy neutrinos and real R with 9i ^ 0, O2 = O3 = and {tun-i , mjv 2 , mjv 3 ) = 
(10 8 ,2 x 10 s , 10 14 ) GeV. Here, tan/3 = 50, Mo = 400 GeV and M 1/2 = 300 GeV. (2a) BR(H X -> rfi). 
Solid, dashed and dashed-dotted lines (the two later undistinguishible here and in all plots) correspond to 
H x — (ho, Ho, Ao) respectively. (2b) BR(fj. — > erf}. The horizontal line is the experimental upper bound. 



The case of hierarchical neutrinos gives clearly larger LFV rates than the degenerate case, 
as can be seen in figs. (J2J), © and ©• However, we will get restrictions on the maximum 
allowed Higgs decay rates coming from the experimental lepton decay bounds. For instance, 
the case of real 9±, that is illustrated in fig. © shows that compatibility with /i — > e7 data 
occurs only in the very narrow deeps at around 6\ = 0, 1.9 and ir. Notice that it is precisely at 
the points 9\ = 0,tt where the BR(Hq, Aq — > r/2) rates reach their maximum values, although 
these are not large, just about 10 -8 . We have checked that for lower tan/3 values, the allowed 
regions in 9± widen and are placed at the same points, but the corresponding maximum values 
of the LFVHD rates get considerably reduced. For the alternative case, not shown here, of real 
02 ^ 0, with 6\ = #3 = we get a similar behaviour of BR{H X — > r/2) with 62 than with 9\, and 
the maximum values of about 10 -8 are now placed at 62 = 0, ir. In contrast, BR(n — > ej) is 
constant with 62 and reach very small values, well below the experimental bound. In particular, 
for tan /3 = 50, M = 400 GeV and M 1/2 = 300 GeV it is lO" 19 . Regarding the dependence with 
#3, not shown here either, a reverse situation is found, where BR(H X — > r/2) is approximately 
constant and, for the heavy Higgs bosons, it is around 10~ 8 . On the contrary, BR(fi — > &y) 
varies but it is always well below the experimental upper bound. In addition, we have checked 
that the BR{t — > /ry) and BR(t — ► cy) rates are within the experimental allowed range in 
all cases. In conclusion, for real R we find that the maximum allowed LFVHD rates are at or 
below 10~ 8 . The case of complex R is certainly more promissing. The examples illustrated in 
figs. (jSJ) and @ are for the most favorable case, among the ones studied here, of complex 62 7^ 



BR(H X -» r/2) Si?(r -> M7 ) 




Figure 3: LFV ratios for hierarchical heavy neutrinos and complex R with 02 7^ 0, 0\ — 83 = 0, and 
(m Nl ,m N2 ,m N3 ) = (10 8 ,2xl0 8 , 10 14 ) GeV. Here, tan/3 = 50, M = 400 GeV, and M 1/2 = 300 GeV. (3a) Depen- 
dence of BR(H X —* Tp) with Arg{9i) for \0n \ = n. Solid, dashed and dashed-dotted lines are for H x = (ho, Ho, Ao) 
respectively. (3b) Same as (3a) but for BR(t — > /ry). The horizontal line is the experimental upper bound. 



with 9\ = #3 = and show that considerably larger BR(H X — > r/2) rates than in the real R 
case are found. Regarding the dependence with 82, we find that for the explored values with 
(\92\, Arg(02)) < (3.5, 1), the Higgs rates grow with both |^| arid Arg(02) and, for the selected 
values of the MSSM-seesaw parameters in fig. (J3J, they reach values up to around 5 x 10~ 5 . We 
have checked that the predicted rates for BR(fi — » cy) are well below the experimental upper 
bound, being nearly constant with 62 and around 10 -19 . Similarly, for the r — > e 7 decay. Notice 
that the smallness of these two decays, in the case under study of 62 7^ 0, is not maintained if our 
hypothesis on #13 = is changed. For instance, for #13 = 5°, which is also allowed by neutrino 
data, we get BR{^jl — * &y) ~ 1.8 x 1CP 8 well above the experimental upper bound. Therefore, 
in this case of complex 62 7^ 0, the relevant lepton decay is r — > ^7 which is illustrated in 
figs. © and Q together with its experimental bound. For the set of parameters chosen in 
fig. ©, we get that the allowed region by r — > fj/y data of the (\02\, Arg(02)) parameter space 
implies a reduction in the Higgs rates, leading to a maximum allowed value of just 5 x 10~ 8 . 
The dependence of the LFV ratios with Mq and M1/2 for hierarchical neutrinos are shown in 
fig. We see clearly the different behaviour of the LFVHD and the lepton decays with these 
parameters, showing the first ones a milder dependence. This implies, that for large enough 
values of Mq or M1/2 or both the BR(t — > /ry) rates get considerably suppresed, due to the 
decoupling of the heavy SUSY particles in the loops, and enter into the allowed region by data, 
whereas the BR(Hq — ► r/2) rates are not much reduced. In fact, we see in figs. (@^) and (jU) 
that for the choice Mq = M x i2 the r decay ratio crosses down the upper experimental bound 
at around Mq = 1100 GeV whereas the Higgs decay ratio is still quite large ~ 4 x 10 -6 in 
the high Mq region, around Mq ~ 2000 GeV. This behaviour is a clear indication that the 
heavy SUSY particles in the loops do not decouple in the LFVHD. Notice also that it can be 
reformulated as non-decoupling in the effective H^Tfi couplings and these in turn can induce 
large contributions to other LFV processes that are mediated by Higgs exchange as, for instance, 
r — > /i/i/x. However, we have checked that for the explored values in this work of Mq, M1/2, 
tan j3, R and that lead to the anounced LFVHD ratios of about 4 x 10~ 6 , the corresponding 
BR(t — > /i/i/i) rates are below the present experimental upper bound. 

In summary, after exploring the dependence of the LFVHD rates with all the involved MSSM- 
seesaw parameters, and by requiring compatibility with data of the correlated predictions for 
fi — ► &y, t — > e7 and r — > /ry decays, we find that BR(Hq, Aq — > r/2) as large as 10~ 5 , for 
hierarchical neutrinos and large MguSY in the TeV range can be reached. These rates are close 
but still below the expected future experimental reach of about 10~ 4 at the LHC and next 
generation linear colliders. 
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Figure 4: Dependence with M (GeV) and M 1/2 (GeV) for (m Nl ,m N2 ,m N3 ) = (10 8 ,2 x 10 8 ,10 14 ) GeV, 6j = 
6 a = and tan/3 = 50. (4a) Bi?(tf -> r/Z) versus M (GeV) for M 1/2 = 300 GeV and 9 2 = Tve 0A \ (4b) Same 
as (4a) but for BR(t -> /ry). (4c) Bi?(# r/l) versus M 1/2 (GeV) for M = 400 GeV and 62 = 7re 8i . (4d) 
Same as (4c) but for BR{t -» ^7). (4e) B.R(#o -» rp) versus M = M 1/2 (GeV) for 6» 2 = vre 8 '. (4f) Same as 
(4e) but for BR(r — > /17). The horizontal lines are the upper experimental bound on BR(t — + ^17). 
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